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Abstract 

The conventional approach to perturbative evolution is illegal because 
the expansion in powers of as of the the DGLAP splitting matrix P(z, Os) 
breaks down at small z. The small- a; data for the proton structure function 
F2{x, Q^) and its charm component (a;, Q^) show that a hard pomeron, 
with intercept close to 1.4, must be added to the familar soft pomeron. 
Conventional perturbative evolution may be applied to the hard-pomeron 
component and provides a striking test of perturbative QCD. 



1 Introduction 

A striking discovery at HERA has been the rapid rise of a'^ ^ with at even 
quite small fixed values of : see figure |l|. If one parametrises the rise as an 
efl[ective power 

then the power A(Q^) is found to be significantly greater than the value 0.08 
that is familiar in purely hadronic collisions The value of X{Q^) has been 
extracted]^ by the HI collaboration from their data and is shown in figure ||, 
from which it can be seen to increase with and reach about 0.4 at the highest 
values that have been measured. 

When W'^ '3> Q^, x is small and W"^ ~ Q^/x. Then the effective-power be- 
haviour (0) corresponds to 

J^2(x,Q2)^/(g2)^-A(Q^) (2) 

When they extracted X{Q^) from their data, to make the plot of figure |^, HI 
assumed that the value of X{Q^) at small x is independent of x at each Q^. 
While the data are compatible with this assumption, they do not require it and 
it does not have theoretical justification. Rather, one should parametrise the 
data with a sum of fixed powers of x, whose relative weight varies with Q^. 
In the two-pomeron model I will describe, the resulting A(Q^) has significant 
variation with x, as is seen in figure O. 
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Figure 1: Data[Q for a'^ 'p — (47r^aEM/Q^)-F2 at various values of Q^, together 
with real-photon data 

Even before the HERA measurements, there were predictions |^ that A((3^) 
would be large at high values of . Such predictions arose from two different 
equations of perturbative QCD: the DGLAP equation and the BFKL equation. 
However, it has since been realised that the predictions are not as clean as 
initially had been hoped. 

The BFKL equation is applicable when x is small, and large but not too 
large. It can be made to give values of A(Q^) close to \, which is approximately 
what experiment finds at high . This led to talk of a second pomeron, the 
hard pomeron, with a X{Q^) that might be calculated from perturbative QCD. 
However, it has been realised that this calculation is almost certainly invalid. 
Apart from concerns that the approximations used to derive the BFKL equa- 
tion do not take sufficient account of energy conservation it is not correct 
to suppose that the equation enables A(Q^) to be calculated from perturba- 
tive QCD aloneQ]. While it is conceivable[|[ ^ that there is no such problem 
with purely-hard processes such as 7*7* collisions at very high Q^, for semi- 
hard collisions such as 7*p the BFKL equation inevitably receives important 
contributions from uncalculable nonperturbative effects. 

The DGLAP equation is derived from perturbative QCD and therefore is valid 
at sufhciently large Q^. It couples the variation of the quark distributions 
in the proton to the gluon distribution. When is much larger than all the 
relevant quark masses, the gluon distribution affects the variation of all the 
quark and antiquark distributions equally, and the variation of the gluon 
distribution itself receives equal contributions from all the quark and antiquark 
distributions. That is, the gluon distribution's evolution is coupled to the singlet 
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quark distribution, which is the sum of the quark and antiquark distributions: 

q{x,Q^) ^u + u + d + d + c + c + s + s + .... (3) 

The number of heavy-quark terms that should be included depends on the 
value of W . So the relevant form of the DGLAP equation introduces the two- 
component quantity 



It reads 




Here, P(z,as((5^)) is a 2 x 2 matrix, called the splitting matrix. The DGLAP 
equation becomes very simple if we introduce the Mellin transforms of u(a;, Q^) 
andP(z,a,(Q2)) 

u{N,Q^)^ [ dxx'^-^u{x,Q^) (6) 
Jo 

and 

P(7V,a,(g2)) = /dzz^P(z,a,(Q2)). (7) 
Jo 

Then 

^u(iv, g2) = p(iv, a,(g2)) u{N, g2). (s) 

In order to apply the DGLAP equation, one should first choose some starting 
value Qq of and fit the parton distributions there to experimental data, as 
functions of x. These functions cannot be calculated from perturbative QCD, 
but the DGLAP equation, together with assumptions about the form of the 
gluon distribution, determines how they change as increases. If[|lO| one 
starts at some fairly small value of with parton distributions that are 
rather flat in x at small x, then an application of the DGLAP equation results 
in distributions that rapidly become steeper as increases, in good agreement 
with experiment ||ll|, However, this application is controversial. It is not safe 
to apply a perturbative equation for choices of Qq as small as 1 GeV^ or less, 
as is sometimes done. Indeed, it is necessary at present to expand the splitting 
matrix P(z,Q;s(g^)) in powers of as and this is unsafe when x is small, as I 
explain below. 



2 Regge approach 

A powerful approach to the small- a; data for F2{x,Q^) is to extend the Regge 
phenomenology that is so successful for soft hadronic processes |l3[ Regge 
theory is not supposed to be a competitor with perturbative QCD, but to co- 
exist with it, and to be applicable when is much greater than all the other 
variables, in particular when it is much greater than Q^. When x is small 
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~ Q^/x, so the condition that W"^ ^ is just that x <^ 1, however large 
may be. 

Regge theory relates high-energy behaviour to singularities in the complex an- 
gular momentum plane ]T^. For deep inelastic scattering, the complex angular 
momentum I is essentially the Mellin-transform of (^. The correspondence 
is 

iV ^ / - 1 (9) 

so that this assumption is equivalent to assuming that the relevant singularities 
in the complex A^-plane are simple poles. The assumption may be not literally 
correct, but it turns out to give an excellent description of the small- a; dataJTsf. 
A pole at = e contributes a power x*^ to the small- cc behaviour of ^2- So 
the soft pomeron contributes x~'^'^ with ei = 0.08. But this is not sufficient 
to describe the rapid rise with 1/x seen in the data at small x and large Q^. 
Another term x~'^° is needed, with eo « 0.4. We call this A^-plane or /-plane 
singularity the hard pomeron. This does not explain what is its dynamical 
origin: maybe it is perturbative QCD though, as I have explained, initial hopes 
that it might be derived from the BFKL equation now seem unlikely to be 
realised. 

So the simplest fit to the small- a; data corresponds to 

F2{x,Q^)^ MQ')^-''- (10) 
1=04 

where the i = term is hard-pomeron exchange and i = 1 is soft-pomeron 
exchange. Regge theory gives no information about the form of the coefficient 
functions fi{Q^), beyond that they are analytic functions. Also, because the 
real-photon cross section is 



47r2a 

g2 ' 



(11) 

Q^=o 



at fixed W the function F2 vanishes linearly with a.t = 0. If we assume 
that each term in ( |lO| ) has this property, then 

MQ') - (Q')'+^* * = 0, 1 (12) 

near = 0. One might hope that perturbative QCD will determine how the 
two coefficient functions behave at high Q^, but as I will explain so far the 
theoretical difficulties allow this only for the hard-pomeron coefficient function 
/o(Q^). Figure ^ shows how /o and /i vary with if we fit (|lO|) to the 
data at each Q^, including values of x up to 0.02. This provides a guide on 
the likely functional forms of the coefficient functions. The fi{Q^) can then be 
parametrised with suitable functions, whose shapes resemble the data points in 
the plots and which include a number of parameters. These parameters, and 



the power eo, may then be varied] 16 so as to give the best fit to all the small- a; 



data for F2{x, Q^) together with the data for a''^. This results in the value 

eo = 0.437. (13) 
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with 



and 



MQ^) = 0.0015 — (14) 

^ ' (l + Q2/9-ll)i+^^« ^ ' 



The fit was made[p^ imposing ei = 0.08 and using only data with x < 0.001 
and therefore < 35 GeV^; see figure ^j. With the addition of a term from 
(/2, 0,2) exchange, and including also powers of (1 — x) in each term to make the 
structure function vanish suitably as x — > 1, the fit agrees [^| remarkably well 
with data at larger x, up to — 5000 GeV^. 

Having concluded that the data for F2 require a hard-pomeron component, it 
is necessary to test this with other data. The hard pomeron is seen clearly 
in the charm component of i^2- This describes events in which a D* particle 
is produced, which are used to extract the contribution F|(a;,(5^) to the com- 
plete F2{x, Q^) from events where the 7* is absorbed by a charmed quark. The 
dataQ for F|(a:,Q^) must be treated with some caution because the experi- 
mentalists have to make a very large extrapolation to compensate for limited 
acceptance. Nevertheless the data have the striking property that, over a wide 
range of Q^, they behave as a fixed power of x: 

F^{x,Q') = UQ')x-''> (16) 

with eo ~ 0.4: see figure |[ It seems, therefore, that F2{x, Q"^) at small x receives 
a contribution from the hard pomeron and that the soft-pomeron contribution to 
it is negligibly small. To a very good approximation the hard pomeron seems to 
be flavour-bhnd[p^ : i^2^(a;, Q^) is close to | the hard-pomeron part of F2{x, Q^). 

For the ranges of x and where they overlap, the two-pomeron fit and con- 
ventional fits based on two-loop unresummed perturbative QCD agree with the 
data equally well. However, they no longer agree when they are extrapolated 
to smaller values of a;, as is seen in figure 0. 



3 Real photons: a crucial question 

Because the data of figure || include a point at extremely low Q^, it should be 
reliable to extrapolate these data to = 0. This extrapolation is the upper 
curve shown in figure || and, at -^/s = 200 GeV, the hard-pomeron component 
contributes about 20 /xb. That is, at this energy, the fits with and without a 
hard-pomeron component differ by about 10%. The errors shown on the data 
in figure ^ are purely statistical; there is an additional systematic error of about 
10% at the lowest Q^. So at present it is not possible to decide whether the 
hard pomeron is present in the a'^P data, though it seems likely. 
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There is also uncertainty with the data for a^'^ from LEP. The cross sections 
are rather sensitive to the Monte Carlo model used for the unfolding of detec- 
tor effects, different Monte Carlos producing different results. In figure ^ the 
resulting uncertainty is contained in the errors on the OPAL data. The L3 
data are shown with the use of two Monte Carlos, the errors corresponding to 
the statistical and systematic errors combined in quadrature. The figure shows 
the prediction obtained by applying factorisation to a^^^ and a^^ data, with no 
hard-pomeron component. The data from L3, particularly the upper set, may 
require such a component, and the figure also shows how adding it in might 
change the fit. 

This is an important question. Is the hard pomeron already present at = 0, 
or is it rather generated by perturbative evolution? If it is there already at 
= it presumably arises because the photon has a pointlike component, as 
if it is present in the pp total cross section it is very small. Maybe it will be 
significant at LHC energies. 



4 Perturbative evolution 

I now show that perturbative evolution governs how the hard pomeron's con- 
tribution to the structure function increases with It is found pt t hat the 
parametrisation of the hard-pomeron coefficient function /o(Q^) of (10) agrees 



very well with what is obtained from DGLAP evolution, over a large range of 
. As yet, we are not able to apply perturbative evolution to the soft pomeron. 

A power contribution f{Q'^)x~'^° to F2{x,Q'^) corresponds to a pole 



in u(7V, Q )■ With four active quark flavours and a flavour-blind hard pomeron, 
fq{Q^) = ^/o(Q^). We find(20|, on taking the residue of the pole at iV = eo on 
each side of the Mellin transform (^) of the DGLAP equation, 

^f(Q2)^p(7V^e,,a^(Q2))f(Q2) (18) 

If we include four flavours of quark and antiquark in the sum in (|^), then at 
= 20 GeV^ the singlet quark distribution xj^filf + Qf) ^ O.OQSa;^*^" at 
sufficiently small x. According to figure ^, the charmed-quark component F2 
of F2 is governed almost entirely by hard-pomeron exchange at small x, even 
at small values of Q^, and, within the experimental errors, its magnitude is 
consistent with the hard pomeron being flavour-blind. According to perturba- 
tive QCD, the charmed quark originates from a gluon in the proton, andpl| 
the two distributions are proportional to each other to a good approximation 
over a wide range of x and Q^. This implies that the gluon distribution also 
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is hard-pomcron dominated. At (5^=20 GeV^ and x = 0.01, a next-to-leading- 
order fit(22[ m to the combined ZEUS and HI data gives xg{x, Q^) ^ 5.7±0.7. 
Other authors ||2^, ^ find much the same value. This is 8 ± 1 times the 
hard-pomeron component of the singlet quark distribution. 

An unresummed perturbation expansion of the splitting matrix P(iV, a^) is not 
valid pO|] for small values of N , but we need P(A^, as) aX, N = 0.437, that is far 
from 0, and so it is reasonable to hope that resummation is not needed. The 
numerical values of the elements of the splitting matrix P(iV, as) are known^ in 
one and two-loop order, so it is straightforward to integrate (|8|). At the energies 
being considered it is necessary to use four flavours throughout as the charm 
contribution is active. The beauty contribution is so small that its omission has 
a negligible effect. 

The result of integrating the differential equation (^8|) for the singlet quark 
distribution is shown in figure where the solid curve is the result of the two- 
loop-order perturbative QCD evolution according to (^8|), and the broken curve 
is the Regge fit to the data I have described. The ratio of the gluon distribution 
to the hard-pomeron component of the singlet quark distribution is taken to be 
8.0 at — 20 GeV^. Figure |lO| also shows how the gluon distribution 

xg{x,Q^)^ fg{Q^)x-'° (19) 

evolves. Provided one chooses A such that as(M^) — 0.116, which is the HERA 
value there is little difference between the leading-order and next-to- 

leading-order results. 

The conventional approach to evolution expands the splitting matrix P(7V, a.,) 
in powers of as- Because an unresummed expansion that needs the splitting 
matrix at small N makes the splitting function larger than it really is, a gluon 
distribution of a given magnitude apparently gives stronger evolution than it 
really should. That is, the conventional approach will tend to under- estimate 
the magnitude of xg(x, Q^) in certain regions of {x, Q^) space. This is verified by 
the results for the evolution of xg{x, Q^) obtained from integrating ([l|). Figure 
|lT| shows the proton's gluon structure function at two values of Q^, according 
to the solution of (|T^), which does not use the splitting matrix at small N ^ and 
compares it with what is extracted from the data by conventional means. 

The agreement between the extraction of the hard-pomeron component of F2{x, Q^) 
from experiment, and its calculated evolution, is a striking success both of 
the hard-pomeron concept and of perturbative QCD. We cannot apply a sim- 
ilar analysis to the soft pomeron, because this would need the splitting ma- 
trix P(iV, tts) at iV = 0.08, which is too small for an expansion in powers 
of as to be meaningful. We do not yet know how to perform the necessary 
resummation p^, |30| , |3l[ |. 



A fuller account of pomeron physics and QCD may be found in a book||3^ to 
be published later this year by Cambridge University Press. 
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Figure 2: The effective power \{Q^) ofdll) extracted from HI datal 
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Figure 3: Data^ for tfie effective power of \/x at = 2.5, 15 and 60 GeV^. 
Tfie liorizontal lines correspond to tlie values plotted in figure |^ and the solid 
lines to the two-pomeron fit. 
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Figure 4: Fits to the coefficient functions fa{Q^) and fi{Q^) of (|ic| ) extracted 
from HI and ZEUS data. The black points are for eg — 0.36 and the white 
points are for eg = 0-5. 
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Figure 5: Regge fit to data for ^2(2;, Q^) for Q2 between 0.045 and 35 GeV^ 
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Figure 6: Data for F^{x,Q'^) from ^ i g to 130 GeV^. The lines are the 
hard-pomeron component of the fit to F2{x, Q^) shown in figure ^, normaUsed 
such that the hard pomeron is flavour-bhnd. 
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Figure 7: Two-pomeron fit to F2{x,Q'^) at various values of (thick curves), 
with two-loop unresummed perturbative-QCD fit|l^ (broken curves) 
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Figure 8: Data for a'^'P . The upper curve is the extrapolation to Q'^ = of the 
curves in figure including also a contribution from (/2, 02 exchange; the lower 
curve omits the hard-pomeron term. 
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Figure 9: Data for a'^^ . The lower curve is the soft-pomeron/reggeon contribu- 
tion; the upper curve has an additional hard-pomeron term. 
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Figure 10: (a) Next-to-leading-order evolution with A = 400 MeV of the hard- 
pomeron coefficient foiQ"^) (solid curve) and the fit to the data (broken 
curve), and evolution of the gluon structure function xg{x, Q^) at a; = 0.01. 
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Figure 11: Gluon structure function xg{x, Q^) at (a) = 20 and (b) 200 GeV^. 
In each case the thick hne is our evolved distribution. In (a) the thin hues are 
the hmits extracted by conventional NLO analysis of HERA data|^, In (b) 
the middle line is@ ^ CTEQ5M and the lower hne is|2|, |2^ MRST20011. 
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